Adaptive hp refinement for spectral elements in numerical relativity

When a numerical simulation has to handle physics problem with wide range of time-dependent length scales, dynamically adaptive discretizations can be the method choice. We present major upgrade relativity code bamps in form fully adaptive, physics-agnostic hp refinement. describe foundations mesh refinement context spectral element methods, precise algorithm used perform bamps, as well several indicator functions drive it. Finally, we test performance, scaling, and accuracy treating 1D 2D example problems, showing clear improvements over static configurations. In particular, consider simple nonlinear wave equation, evolution real scalar field minimally coupled gravity, gravitational waves.